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As a continuation of the work in we discuss the Casimir effect for a massless bulk scalar field in 
a 4D toy model of a 6D warped flux compactification model, to stabilize the volume modulus. The 
one-loop effective potential for the volume modulus has a form similar to the Coleman-Weinberg 
potential. The stability of the volume modulus against quantum corrections is related to an appro- 
"jC^ ' priate heat kernel coefficient. However, to make any physical predictions after volume stabilization, 

knowledge of the derivative of the zeta function, £'(0) (in a conformally related spacetime) is also 
. required. By adding up the exact mass spectrum using zeta function regularization, we present a 

■ revised analysis of the effective potential. Finally, we discuss some physical implications, especially 

concerning the degree of the hierarchy between the fundamental energy scales on the branes. For a 
larger degree of warping our new results are very similar to the ones given in Ref [lj and imply a 
larger hierarchy. In the non-warped (rugby-ball) limit the ratio tends to converge to the same value, 
independently of the bulk dilaton coupling. 
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Recent studies on braneworld models with extra dimensions that are compactified by a magnetic flux and bounded by 
codimension 2 branes (conical singularities) have revealed various fundamental properties, i.e., the regularization 
of branes (conical singularities) /linearized gravity 0, H, @| and stability against classical perturbations @, [1] . There 
are still a lot of unsolved questions which should be clarified, especially concerning non-linear gravity and cosmology. 
Furthermore, careful analysis of these simpler models may give important physical insights on self-gravitating branes 
in various flux compactifications in string theory Q . 

In this paper, we focus on models of braneworld based on a 6D supergravity though we will work with its 
4D counterpart [l(|, because of the lack of a formulation of the heat kernel coefficients for conical singularities in 
6D. In these supergravity inspired models, the size of the compactified internal-space is not fixed classically and may 
behave as a volume modulus in a 4D effective theory. To fix this size of the volume (the volume modulus), we should 
therefore discuss other additional mechanisms. The authors of Ref. [l[ focused on Casimir corrections for the case of 
perturbations of a massless, minimally coupled bulk scalar field. In this article, we present a revised analysis of the 
previous calculations, by performing a strict mode summation of the exact mass spectrum, which is now wholly taken 
into account. 

We give an exact analysis of the one-loop effective action in a 4D alternative model. The one-loop effective potential 
for the volume modulus can be written in a form which is somewhat similar to the Coleman-Weinberg potential: 

Vi,eS = , U-l) 

P+ 

where p+ characterizes the volume modulus and A± and B± are functions of the model parameters and the shape 
modulus (the degree of the warping of the bulk, which is completely determined if one fixes the brane tensions). 
Stability itself is determined by the sign of B4, which is closely related to the relevant heat kernel coefficients. 
However, to discern phenomenological effects on the brane, i.e., the effective mass of the modulus and the degree of 
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the hierarchy between fundamental energy scales on the brane, we need to know the value of A4 as well. A4 is not just 
related to the heat kernel, although partly related with it; we also need to evaluate C'(0). Thus, an accurate evaluation 
of A4 is crucial to make physical predictions, e.g., for the hierarchy or cosmological constant (CC) problems. 

In the work 1], it has been shown how to divide A4 into two parts, by introducing a continuous conformal trans- 
formation, going from the unwarped frame to the original warped geometry The first part is associated with a one 
parameter family of conformal transformations and is called the cocycle function. This term can be obtained from 
heat kernel analysis. The second part consists of calculating the derivative of the zeta function in the unwarped frame, 
where to evaluate this piece we used the WKB approximation. In this paper we shall present an exact analysis of the 
mass spectrum for these Kaluza-Klein like modes in the unwarped frame. 1 

The structure of the article is as follows: In the next section we briefly discuss our 4D analogue warped flux 
compactification model and review the tools to analyze the one-loop Casimir effect. We give the exact mass spectrum 
in the unwarped rugby-ball frame. In Sec. Ill we carry out the zeta function regularization in the unwarped frame. 
In Sec. IV, we give physical implications of our result, e.g., for the hierarchy problem by making comparisons with 
the results in PJ. We finish this article after giving a summary and some discussions. 



II. 4D WARPED FLUX COMPACTIFICATION MODEL 

Our main interest is the Casimir effect in the warped flux compactification model in 6D supergravity Q : 

5 6 = Ml J d 6 x^ (~R - \d AV d\ - \e-*F AB F AB - 2g 2 e^ , (2.1) 

where F A b is the field strength of the electromagnetic field, ip is the dilaton, and go is the bulk dilaton coupling. 
Hereafter, we set M| = 1 and if needed, we put it back explicitly. In order to evaluate the Casimir energy in this 
model, we need to know the heat kernel coefficients on the conical singularities. However, as mentioned in there 
is no mathematical formulation of these coefficients the authors are aware of. Thus, in this paper, we focus on its 4D 
counterpart theory: 

S 4 = M 2 f d*xV=g (r - \d AV d A ^ - \e^F AB F AB - 4g 2 eA , (2.2) 



2 

which has a series of warped flux compactification solutions [l(]| 

ds 2 = h{p)d9 2 + ^ + (2p)(-dr 2 + dx 2 2 ) , h(p) = ^( P+ - p) (p - p.) , 

Hp) p 

<p(p) = -]n(2p) , F Bp = -%. (2.3) 

P 

We also set M4 = 1 unless it is needed. Then, the branes correspond to strings, lying at the positions given by 
the horizon-like condition h(p) = 0. It is useful to define the coordinate z = {{p+ — p)/(p — P-)) 1 ^ 2 , to resolve the 
structure of spacetime in the case that p+ = p_. For z — > 0,oo, there are cones whose deficit angles are obtained 
from the relation S± := 2ir — g 2 {p+ — p-)A8/(p±). 

The moduli in the effective theory on the branes are characterized by p±, or equivalently by p+ and r, where the 
string tensions are related to the conical deficits by a± — M\b~±. We stress that these relations are only valid for 
sufficiently small brane tensions, in comparison with the bulk scale M 2 . The coordinate 9 is defined as < 8 < A8 
and A9 is given by 

M(r,6 + ) = ^-^±- (2.4) 
g 2 (l - r) 

where r :— p-/p + . Once the brane tensions, er+ and er_ are fixed, then r is also fixed and so we now regard as free 
parameters r and 5 + , along with the dilaton bulk coupling g. The remaining degree of freedom used to determine the 
bulk geometry is the absolute size of the bulk, i.e., p+, which should be fixed by additional mechanisms. 



1 These modes are not the standard KK modes, because of the presence of conical singularities at the poles of the two-sphere on the 
internal dimensions. Given the similarity to a rugby ball we shall also call this unwarped frame the rugby-ball frame. 
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We discuss the basic properties of a massless, minimally coupled scalar field on this background: 2 



Sscala 



Now consider a continuous conformal transformation of the metric, parameterized by e; 

dsl^e^-^dsi, w=iln(2p), 



(2.5) 



(2.6) 



where for e = 1 we have the original metric, which we shall define as A4. e = A4, and e = is the conformal frame. 
The classical action changes as 



^scalar = ~ J d A Xyfg(j)A 4 (j) = — - J d i Xy/§4> (A4 + E 4 (e) 



(2.7) 



where 



E 4 (e) = -(e-l) 2 g ab W a ujW b uj + {e-l)A 4 lnuj 

= ( 1 V ff 2 (l - e)(p+ P-){p+(2 + (1 - e)z 2 ) + p_z 2 (-l + e 2z 2 )} 
\2p) (P++P-Z 2 ) 2 

We now derive the mass spectrum in the unwarped frame: 

(a 4 + £(0)) 4>x = -\ 2 4> x , 
We shall decompose the mass eigenfunction as 

The equation of motion of equation (|2.9p has a series of exact solutions : 

^2 \ — n/2K / ^ \ nr /2^ M ^ 



(2. 



(2.9) 



(2.10) 



$ A (r) = 



1 + rz 2 
1 + z 2 



A 



1 + z 2 



D 



2 \ u/2k 



1 + z 2 



1 + Z' 



1 + z 2 

— nr /2k 



i F,(l-,--(l- r ,,,--,l-,.).l--; r 



(It, , IV, , < (/ \ 

1 - v + ^ 1 - r ^ + ^ 1 - r )' 1+ K ; T+^) 



(2-11) 



where k := 1 



" = 2 1 + V J 



A 2 - k 2 n 2 
+ ^(l-r) 2 
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(2.12) 



Here we assume that the mode functions are regular on both conical boundaries. In order to ensure such a condition 
at z = we must have A = for n > (and B — for n < 0). Thus, we arrive at the following mass spectrum 



2\n 



4n 2 r n 



A 2 =fc 2 + 5 2 [4m(m + l) + -^(2m+l)(l + r) + — ], m = 0,1,2,- ••. 



(2.13) 



Note that this mass spectrum is very similar to the case of the gauge field perturbations [l2j ■ (See also @ for exact 
solutions of massless tensor perturbations). 3 



2 A related model including a self-interaction term was discussed in 

3 Our method to determine the mass spectrum is essentially based on the same arguments as that in [l2| . They impose two independent 
conditions, namely normalizability and regularity (denoted as "hermitian conditions" in [I2I0 . However, if we were to only impose the 
condition of normalizability then there may still be other KK-type modes, though they are somewhat out of the scope of this paper. 
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From the above analysis, we can compare the results of the WKB approximation with this exact one. Note that the 
dependence on the winding number n appears only in the form of \n\. This fact just means that the internal space is 
axisymmetric. This mode sum will be the starting point for our analysis of the zeta function. 

The coefficients in the one- loop effective potential, Eq. (jl.ip . can be written in the form 



d 2 xA 4 (r, 5+) 
d 2 xB 4 (r, S + ) 



d 2 x^ 5 " 



P+ 



de a4 (/ = iln(^))-k'(0 ! A 4 , e=0 ), 
2 P+ 2 



2 _B 4 (r,S + ) /vn a a 1 ft i\ 

d x = -C(0, A 4 . e=0 ) = 77 °4 / = 1 , 

P+ 2 2 



where the zeta function £(s) is defined for the given mass spectrum in the unwarped frame Eq. (|2. 13|) : 



(2.14) 



(2.15) 



The coefficients other than — £'(0)/2 in the conformal frame were derived in Ref. [l|, where it was shown that the 
term A 4 (r, <5+) is given by the equation 

Mr, S+) = (Mr, S + )) cocyclc + (Mr, M)_ ped > (2-16) 

where the first term is the cocycle correction. In the next section, we will discuss the evaluation of the second term, 
which is — £'(0)/2 in the unwarped, factorizable frame: 



d 2 i 



1 



= --C'(0,A 4 , £=0 ). 

unwarped Z 



{Mr,5 + )) 

The cocycle function can be derived from heat kernel analysis |l|: 4 

)tt 2 J J \2 \p + JJ (l + z 2 )(l 



(2.17) 



^4, cocycle <^+) — 



14407 



+ rz 2 ) (1 + rz 2 ) 4 



9 2 S + 2-^t 
144tt 2tt I _ £+ 

144tt r 1 _ i± 



1 1 + (1 _ £t)2 n 
21n(2)-l + r-g l _ S J*J (- h(2)(l - 3r) 



1 - r 



1 lr 2 + fl-^i) 2 /l 
2 ln(2r) - 1 + - - - 1 ( - ln(2r)(l 

V ' r 5 (l_it)2 V2 V A 



1 - 



(2.18) 



where ^(p+, p_, e, z) can be found in Appendix B of [l| and is just a function of r and z. The independence of $ on 
e is discussed in PJ. The coefficient B 4 (r, 5+) is given by 



where 



and 



B 4 {r, ^ 



(Mr,S + )) hul] 



B 4 (r, <$ 4 



g 2 (27r-^+) 
1440tt 2 



bulk 



B 4 (r,6 + ) 



branc 



dz 



*(l,r,0,z) 



(l + z 2 )(l + rz 2 ) (1 + rz 2 ) 4 ' 



■^-^^r4-i 1 + ( 1 -fe) a (i- 3 r) 
288tt 2tt i _ *£ L 5 (i _ |^) 2 

<5 



(2.19) 



(2.20) 



2ti 



2 5 

9 

"288tt 



§i - (1 - r) r + (1 - §i) r 1 r 2 + (1 - |±) 2 , 3 N 

L 5 (i-§+)2 { r 



(2.21) 



which is also derived from heat kernel analysis. 



4 See also |13H . for UV effects and heat kernel coefficients in higher dimensions. 
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III. ZETA FUNCTION REGULARIZATION IN THE UNWARPED CONFORMAL FRAME 

Given the eigenvalues we found in the previous section, we shall now derive an expression for the corresponding 
spectral zeta function and related quantities, such as the effective action. In Ref [l| the authors used the density of 
states method because they had not found an exact mode spectrum. 

The first step on our journey will be to use the Mellin transform: 



C(s) = — r-r ( d 2 x /d 2 fcV f i s_1 exp{-[fc 2 + g 2 (a(m + a) 2 + b(m + a)\n\ + cn 2 + q)]t} . 



(3.1) 



where we have just rearranged the eigenvalue X mn k , see equation (|2.13|) , such that it takes the form of a two dimensional 
Epstein zeta function, and 

4fl+r1 4r 
a = 4, 1 j , c=^, q=-l, a = 1/2. (3.2) 

K K 

Firstly, we perform the fc-integration after interchanging the order of integrations: 

/ d 2 x f d 2 kJ2\- 2 ; k = -|- / d 2 x f d 2 A;V I dtt s - 1 exp{-[k 2 + g 2 (a(m + a) 2 + b(m + a)\n\+cn 2 + q)}t} 
J J m,n l( - S ' J J m ,n Jo 

f 27T f°° 

= d 2 xy2—— dtt s ~ 2 c?vp{-g 2 [a(m + a) 2 +b(ni + a)\n\+cn 2 +q]t} 
J m,n T ( S > Jo 

= J d 2 x ^j^- Y} a ( m + a ? + b(m + a)\n\ + cn 2 + g] 1 ~ s (3.3) 



where in the third step we have redefined the terms such that 

A/ , bn 1 (l + r)|d 

a(n) = a H = - + — . 

V ' 2a 2 2k 

2 (bn) 2 n 2 

4a k 2 



( 7(n)=g + cn 2 -^ = -^(l-r) 2 -l. (3.4) 



This form shows that the zeta function is essentially some kind of two-dimensional Epstein zeta function as we show 
below. In the following discussion, we focus on the non-axisymmetric modes n / 0. For the contribution of the 
axisymmetric modes, see Appendix A. 

A. Extended binomial expansion 

Frequently, the form of the two-dimensional Epstein zeta function allows one to perform its summation in an elegant 
way which involves the Chowla-Selberg expansion formula or more frequently (as it would correspond to the case here) 
a generalization thereof, see [T3|. However, some conditions must be fulfilled, in particular, the quadratic form must 
be positive definite and the constant q term should be also non-negative. But this is not here the case, and the fact 
that q < does not allow for such a beautiful analysis. Henceforth, we shall apply in what follows what we have 
called the extended binomial expansion approach (see also flU). n a- m ely, we will introduce an extra summation via 

oc oo 

\a(m + a) 2 + b(m + a)n + cn 2 + q] 

m— n— 1 

OC OO OC OO OC -p / . \ 

= EE Hm + o) 2 + q] - s+1 = E E E r(s-iw [a{m + &)2] ' (3 - 5) 

m =0n=l m =0n=l j=0 

where the validity of the binomial expansion is defined for 

' < 1 , (3.6) 



a(m + a)' 2 

which is indeed satisfied for all possible values of m and n in our model. 
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The crucial step is to interchange now the summation over j with the m, n summations. The absolute convergence 
of the binomial expansion in j guarantees such an interchange. Whence, we obtain that 

E E H m + ") 2 + 6 ( m + Q ) n + C1? + «] ~ S+1 = E E a 1 - 8 ^ p fs _iV! V Cg(2« + 2j - 2, a) , (3.7) 

m=0n=l j=On=l 

where we have summed over n to give a standard relation in terms of Hurwitz zeta functions. Substitution of the 
above relation into Eq. (|3.3p leads to 

,2(1-) f„4^„, x9 /„ „ „1 l+r 

-n 



cwu = ^/^E^^E&a-o'+i) Cff (2 S+ 2i-2,i 



j"=0 n=l 



2k 



(3.8) 



with 



2-w-Pr( s + j-i) 
G(j ' s) = ■ (3J) 

Note also that the presence of \n\ in our mode sum, with — oo < n < oo, makes the zeta function a non-typical one. 
In our approach, we have just multiplied by two and taken the sum from n = 1 . . . oo, the n — mode having been 
treated independently. 

The above expression does not converge unless Re s > 3 and thus, we must first analytically continue the zeta 
function in equation (|3.8[) in the usual way by first subtracting off the terms responsible for the divergences in 
C(s). After this we are able to take the derivative at s — 0, being careful to add back the terms we subtracted off 
before the analytic continuation. 



B. Analytic continuation of the £ function 

In order to isolate the divergent behavior in the zeta function, due to small s, we need to make an asymptotic 
expansion of Eq. (|3.8|) . the details of which are given in Appendix A. From this we can now analytically continue £(s) 
and take the derivative at s = 0. Subtracting off the divergent terms from the zeta function leads to a function 



2(1 — s) « oo oo 2 -j -1 _i_ 

P(a) = 9 —^- /^E G ^ s )E[(^( 1 - r ) 2 + 1 ) Cff(2 S + 2j-2,- + -^n)-F(n,j; S )], (3.10) 

3=0 n=l 



where P(s) is now a regular function at s — and F(n,j; s) is defined in Appendix A. Performing the sum over j and 
adding back the subtracted terms leads to the analytically continued zeta function 

oo 

COO ^ =P(s)+Y,G(j,s)A(j >S ) (3.11) 

where, given that the subtracted terms are simple functions of n, they can be expressed as a combination of ordinary 
Riemann zeta functions: 

A ,,) . <L— J (1 <_y + > J+ „ ^ + 2 . _ [gfaP. - 3) + ^fi> fa(2« - 1) + »)fa(2 a + 1) 

(3.12) 

the functions ifo, u>2 and W4 being defined in Appendix A. Note, however, that we must leave the sum over n in P(s) 
as it is. The j = and j = 1 terms must be considered separately because of the gamma functions in G(j, s) (see 
Eq. (|3.9p ). which leads to compensations of poles and zeros of the Riemann zeta functions in Eq. (|3. 12[) . Taking all 
of this into account, the analytic continuation of the derivative of P(s) with respect to s can be found without any 
further problems. 

Then, after analytic continuation to s — > 0, we obtain P'(0) as given in Appendix B. In Fig. 1 we have plotted the 
function P'(0) as a function of r. As can be seen, the summation over j converges fairly rapidly and the result for 
Jmax = 10 already coincides well with that for j max = 50. In the following calculations, we shall use j max = 50 as 
a conservative choice. Note that in all of our plots the n-summation converges fairly quickly due to the asymptotic 
expansion and the first few terms would suffice; however, we will always set n max = 50 as a conservative choice. 
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FIG. 1: Plots for P'(0), Eq. flB2j, as a function of r are 
shown for S = 0.01. The red and blue curves correspond 
to the truncation of the j-summation at j ma x = 10 and 
jmax = 50, respectively. 
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FIG. 2: Plots of the derivative (d/ds)(J2j GA)(s)\ s ^ , 
Eq. (|3.13|l . as a function of r are shown for S — 0.01. The 
red and blue curves correspond to the truncation of the 
j-summation at j max = 10 and jmax = 50, respectively. 



The analytic continuation of the counterterms Y]j G(j, s)A(j, s) is discussed in Appendix B, where after considering 
the j = 0, j = 1 and j > 2 modes separately we find (see App. B) 



d_ 

ds 



oo 



3=0 



ds 



G(0,s)A(0,, 



d 
ds 



G( M)A ( M))| s=o + E|(^, S )A ( i !S ))| s=o( 3.13) 



J=2 



In Fig. 2 we plot Eq. (|3 . 13|) as a function of r. The j-summation truncated at j max = 10 already exhibits good 
convergence and coincides well with that for j max = 50. In the following calculations, we use j max = 50 as a 
conservative choice. 

We finish this section by noting that the above, quite involved mathematics have yielded in the end very precise 
numerical results. We will now take advantage of them by considering their implication — e.g., those of the Casimir 
corrections from bulk fields — for the phenomenological predictions that were mentioned in the first section. 



IV. IMPLICATIONS FOR THE HIERARCHY PROBLEM 

We can now consider the physical implications on the brane in the original 6D theory given byEq. (J27TJ). InRef (3, 
several phenomenological applications have been discussed, i.e., to the hierarchy problem and the vacuum (Casimir) 
energy density (i.e., the effective 4D cosmological constant) on the branes. With our new results at hand, we give 
a revised version of them, especially focusing on the hierarchy problem. After analysing the hierarchy problem, we 
briefly discuss the case of the vacuum energy density. 

In the original six-dimensional model, the effective four-dimensional Planck scale is 

M P 2 ^ P+(2 V' +) M 6 4 . (4.D 
9a 

If we assume a brane localized field whose bare mass is given by m 2 on either brane at p± then the observed mass 
scales are m 2 + = m 2 and m 2 _ = r 2 m 2 . Thus, the mass ratio between the field and the effective Planck mass is given 
by 

m± _ [l 2 m 2 __3o_ -(2A 6 +B 6 )/(2B e ) u r,s 

A/ 2 - Mi 2tt-5 + ' 1 ' 

where A$ and B 6 correspond to A4 and B4 in the 6D model, respectively. Assuming the factor (/j,to/M|) 2 takes an 
optimal value of 0(1) at the unification of the fundamental scales, then the mass ratio becomes 

n 2 

gp C -(2A 6 +B 6 )/(2B 6 ) ^ / 4 3 n 

2tt — 5+ 
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FIG. 3: Numerical plots of log 10 R(r, 8+ = 0.01) as a FIG. 4: Numerical plots of log 10 7?(r, <5+ = 0.01), via the 

function of r are shown for g = 0.5, 5, 50 (the red, green WKB approximation, as a function of r are shown for 

and blue curves, respectively). We set j max = n max = 50 9 = 0.5, 5, 50 (red, green and blue curves, respectively), 
as a conservative choice. 



where we have used the value of p+, at stabilization. The corresponding quantity in the 4D toy model is given by 

R( r> 6+ ) = 9 \ e -(A i+Bi )/Bt (4 4) 

ZTT — + 

In Fig. 3 we have plotted log 10 [R(r, 0,01)1 as a function of r. The resulting ratios between the scalar and gravi- 
tational energies lead to the same conclusions as in the analysis given in Ref. [ij , in which the derivative of the zeta 
function was evaluated by the WKB method (See Fig. 4). Namely, wc obtain larger mass hierarchies for smaller r 
and smaller g; in detail, however, the form of R is not the same. The most important difference is the behavior at 
r = 1: in the exact spectrum R does not depend on the coupling g (whereas in the WKB approximation it scales 
as g 2 ). Numerical plots also show that this feature is almost independent of the value of the deficit angle S + . Thus, 
the expected mass hierarchy in the rugby-ball frame is independent of the value of the dilaton coupling g, which is 
actually the bulk cosmological constant, because the dilaton dynamics are absent in the conformal (rugby-ball) frame. 

Finally, we add a few comments on the Casimir energy density realized on the branes, i.e., the effective 4D cos- 
mological constant. As shown in Eq. (56) (also see (57)), the expression for the vacuum energy density has a very 
similar form to i?, since B^r, S+) a j 2 . In fact, after setting /i = M4 as an optimal choice, numerical plots of Eq. 
(57) in Ref. [l[ as a function of r for different dilaton coupling g again show convergence to the same point in the 
r —> 1 limit, though here we do not present these results explicitly. For smaller r, the same behavior as in Fig. 5 
in [l| is obtained and implies a tiny amount of effective cosmological constant on the branes, in comparison with 
the gravitational energy scale Mf. The only problem is that the sign of the energy density becomes negative, as is 
expected from the form of the effective potential Eq. (jl.ip . This fact requires an additional uplifting mechanism for 
the effective potential. 



V. CONCLUSION 



In this article, we have given a precise analysis of the mode summation in deriving the one-loop effective potential 
for the volume modulus in a 4D analogue model of a warped flux compactification based on 6D supergravity. In 
[l[ the mass spectrum was found using the leading order WKB approximation. Strictly speaking, this is only an 
approximation and an exact one-loop analysis of the mass spectrum has been performed in the present paper, by 
employing zeta function regularization techniques. We would also like to stress that the unwarped frame is not the 
standard Kaluza-Klein spacetime, because of the presence of conical singularities. 5 

In this paper we used our results to look at the hierarchy problem in warped flux compactification models. The 
qualitative features in Fig. 3 are similar to those of Fig. 4: we expect a larger mass hierarchy for smaller r (where 



5 A beautiful discussion of zeta function regularization in standard KK theories can be found in [15| . 
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FIG. 5: Numerical plot of the integrand of £(0) and 2Bi 
as functions of g are shown for r = 1 and S+ = 0.01. The 
red and blue curves correspond to the cases of ("(0) and 



2B4, respectively. We set j n 



= 50. 



FIG. 6: Numerical plot of the integrand of £ (0) and 2B4 
as functions of r are shown for g = 10 and 5+ = 0.01. 
The red and blue curves correspond to the cases of £(0) 
and 2B4, respectively. We set j max = n max = 50. 



r is the degree of warping) and for larger bulk dilaton coupling g. Indeed, for r <C 1, the contribution from the 
cocycle function becomes important and gives rise to a large mass hierarchy on the brane, though at the same time 
the effective mass of the modulus may become lighter (and hence more unstable to KK perturbations). However, 
quantitatively there is a significant difference for r < 1, where quite surprisingly all the curves for different values of 
g in Fig. 3 converge to the same point at r = 1. In Fig. 4, which is for the WKB approximation, the behavior of R 
is proportional to g 2 , even at r = 1. This may well represent a non-trivial feature of the Casimir effect in the rugby 
ball limit. Similar results are obtained in the case of the vacuum (Casimir) energy density on the branes, namely 
convergence to the same value in the r — > 1 limit, independently of g. For smaller values of r, we can expect a tiny 
amount of the effective cosmological constant on the branes, in comparison with the gravitational energy scale. The 
negativity of the energy density may require slight modifications of the background model. 

As a consistency check of our results, in Figs. 5 and 6 we have plotted the integrands of C(0) and B^r, <5 + ), which 
are related to a^{f — 1) given in Eq. (|2. 14[> . as functions of g and r, respectively, for a fixed value of the deficit 
angle, <5+. Note that both quantities are practically insensitive to 5 + . In the unwarped frame, due to the properties 
of the heat kernel coefficients (see e.g., [la]), the equation £(0) = 04 (/ = 1) should be satisfied. Although they do not 
agree exactly, they do exhibit a similar behavior, namely increasing for larger g and for smaller r. We thus conclude 
that our results are physically reliable. Nevertheless, there are still slight differences between these values and it is 
worthwhile to consider the possible origin of this difference. 

In [l[ the authors confirmed that B4 (and thus 04(7 = 1)) is independent of the value of e which characterizes 
the conformal transformation given in Eq. (|2.6[) and therefore, we believe that B4 is correct. 6 We have also carefully 
checked the convergence of the mode summations. Thus, it is the view of the authors that a possible origin for the 
disagreement could be in our determination of the exact mode spectrum. By imposing vanishing conditions for the 
mode functions and their derivatives at the poles (branes), we were able to derive the mass spectrum in the rug by-ball 
frame. Our method to determine the mass spectrum is essentially based on the same arguments given in [121 ] , where 
two conditions, normalizability and regularity at the poles were imposed (more precisely, it was demanded that the 
wave functions should be Hermitian). However, if only normalizability were imposed, this would allow for logarithmic 
divergences at the poles and this may well lead to additional modes in the eigenvalue spectrum, which could then 
reflect in the differences between Figs. 5 and 6. Regardless of this, based on the qualitative behavior of B4 and C(0), 
the existence of such modes does not seem likely to significantly affect the qualitative behavior of our results. 

In closing we would also like to stress that although the results presented here are strictly for a toy 4D warped flux 
compactification model, the methods introduced can be generalized straightforwardly to six dimensions. However, as 
discussed in [l|, to start with the relevant heat kernel coefficient on the cone in six dimensions should be found. In 
addition, in the current investigation we have focused on the simplest case: a massless, minimally coupled scalar field. 
It should be also possible to extend our analysis to the case of a bulk scalar field with self interactions and other fields 
in the multiplets appearing in the supergravity model, which our model is based on. We hope to report on the results 
of these extended analyses in the near future. 



If any terms in B4 were incorrect then conformal invariance would be ruined. 
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APPENDIX A: ASYMPTOTIC EXPANSION OF THE ZETA FUNCTION 

In order to isolate the divergent behavior in the zeta function, due to small s, we need to expand Eq. (|3.8[) . Towards 
this end we shall make use of the asymptotic expansion of the Hurwitz zeta function [17j ] : 

1 1 + r N 1 (1 (1 + r)n\-( 2s +2J-3) \ ,\ (\ + r )n\ -(2*+2j-3)-i 

Ch(2s + 2j-2,- + — — n) - ' ■ ' 1 » 



2 2k ' 2s + 2j - 3 V2 2k J 2 V 2 2k 

J fc=2 V 7 

where the -Bfc are Bernoulli numbers. Furthermore, in the limit of large n and fixed j, the asymptotic expansion of 
£h is found to be 



1 1 + r 



j 



( H (2s + 2j-2^ + ^-—n) ^(l- r ) 2 + l ~ F(n,j;s), (A2) 

A ZiK L K 

where 

F (nJ;s) := — _ 45(2j + 2j ._ 3) ( -n- + ) . (A:!) 

and we have defined 

w (j,s) = 360(1 + r) 4 , 



360j(l + r) 4 

v l -r) 2 

180(j - l)j(l + r) 4 60j(l + r) 2 (-3 + 2j + 2s)(-2 + 2s + 2j) 



W2 (J' S ) = ""n^ 32 ' - 6 °( 1 + r ) 2 (-3 + 2s + 2j)(-2 + 2s + 2j) , 



(1-r) 4 (1-r) 2 
+ 7(-3 + 2s + 2j)(-2 + 2s + 2j)(-l + 2s + 2j)(2s + 2j) . (A4) 

Note the form of these terms, especially 11)4, will be important for the analytic continuation of the zeta-function. 

APPENDIX B: ANALYTIC CONTINUATION 

In this appendix we discuss the analytic continuation of the various functions discussed in the main text. The 
expression for the analytic continuation of £(s) to s — ► is given by 

2 r 00 1 1 _i_ 1 i 2 

C(0) = 9 - J d 2 ^[-4(c ff (-2,- + ^n)-F(n,0,0)) + (cH(0,- + ^n)(^(l-r) 2 + l)-F(n,l,0)) 



n=l 

«5 2 /\ 2 r/(l + 4 (1+r) 2 „) (l-r) 2 /3(l + r) 4 3(1 + r) 4 2 



7207r(l + r) J ™ ~W k 4 k 2 J (1 + r) 2 V k 4 + (l + r) 2K U 

^2, (l-r) 2 ^ / (1 + r) 4 (1 + r) 2 \1 
+ g (2J - 3 )(l' + rfl ( m JT^W - 6 °?T^F (2J - 3> + MW - 3 " 2J - ")] • <B1 » 

Note that the axisymmetric modes do not contribute to C(0)- 
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We next discuss the analytic continuation that is required to derive the expression for C'(O). The analytic continu- 
ation of P'(s) to s — > is given by 

P'(0) = M{( 2 ln(2fl)-l)f;[c ff (-2^ + ^n)-F[n,0;0]' 



2k 



ln(2 5 ) 



1 °° 1 1 _L 2 

+ 4 E [2C*'(0. 2 + - r ? + X ) - ^ s 

^ 2" 2 ' ^ / fn 1 1 + r ,/n 2 , 2 



^Ki;0])} 



(B2) 



Note, we must also pay attention to the counterterms. The product G(j, s)A(j, s) can be written as 

1 T(s + .]-!) g 2 ^K 1+2s {l - rfi 



G(j,s)A(j,s) = 



— r / d 2 x 

s + 2j - 3) J 

( woO,s) j \-w 2 {j,s) - \-W4{j,s)Cr{2s + 1) , 



2 4 tt r(s)j! ( 1 + r )i+2 S+2j 4 5 ( 2 



(B3) 



where, as before, we must consider the j = 0, j = 1 and j > 2 modes separately. The contribution from the j = 
mode is the simplest one, and analytic continuation to s = results in 



1( G (0, S )A(0, S ))| 



1 1 g 2 n 



+ o 2 4 tt 135 (1 + r 
,C*(-3) 



(B4) 



(wo (0,0) 



K 

2Cfl'(-3) 



^2(0,0) ^2 + w 4 (0,s)Cfl(2s + 1) 



s^O 



- Wo (o,o) ^,; 4 " y +W (o,o)^^ + 2 W2 (o,o) ^ 1} 



ft 



+ d s ( W4 (o,s)G?(2s + i))| s } 

- ^-T^T /^^ + 21nf^))f3ii±^ + 30il±^-42) 
2 4 tt 135 (1 + r) J IV3 Vg(l + r)//V k 4 k 2 7 

(1 + r) 2 720(1 + r) 4 720(1 + r) 2 
~ 50 — ~i — + Zl & (~ 3 > Z2 Cfl (-1) + 14(11 - 67) \ , 

where we have used the fact that 

w 4 (0, s)( R {2s + 1) = -42 + 14(11 - 67)s + 0(s 2 ) . (B5) 

Similar steps follow for the j = 1 mode. By setting j = 1 in G(j, s)A(j, s), differentiating with respect to s, and 
then taking the analytic continuation as s — > 0, we obtain 



^(G(M)A(M) 



1 1 g 2 n{l-rf 



,0 2 4 tt 45 (1 + rf 

,(l + r) 4 



/^{( 2+2 '°fe)) 

60(1 + r) 2 ^ 



-10 



»(l-r) 
(1 + r) 2 720(1 + r) 4 , 



(1-r) 2 



+ 



Cfl'(-3) + 



720(1 + r) 4 , 



-14 



K" 1 k 2 (1 — r) 2 

120(1 + r) 2 / 120(1 + r) 2 



Ch'(-i) 



(1-r) 2 



+ -'(-- + W)}' 



(B6) 
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where we have used the expansion around s = 0: 
w 4 (l, S )( R (2s + 1) = -14 + 
For the modes with j > 2, we obtain 

E^(G(i, S )A(i, S )) 

^ oo 
94^E 



-14- 



120(1 + r) 



^H 8 + 



120(1 + r) 2 
(1-r) 2 



s + 0(s 2 ). (B7) 



2% ^ 45j(j - l)(2j - 3) (1 + rY+V 

3 — z 



\d l + r)J\ V 



'w 4 (j,s)( R (2s + 1) 



x {[^'-^-y-a '-V fl( i + r);JV r (s) 

+ w (j, 0)^^ + W2 (. ? , 0)^^ + \w 4 {j, 0) + ]- (3 1W4 (j, 0) + d s w 4 (j, s) ) } 



2 4 ^ ^ 45j(j - l)(2j - 3) (1 + r)!+^ 
{[«(i-l)--^ + 21n( 



^(1 + 



x 2j(i-l)(7(-3 + 2j)(-l + 2j) 



30(-3 + 2j)(l + r) 2 45(1 + r) 



+ 



+ 



3(1 + r) 4 5(1 + r) 2 [6 + 4 ? 2 - 8 j 2+r ^ 1+ ^ 



+ ■ 



(1-r) 2 

(1-r) 2 



+ 



(1-r) 4 



K 

2 



(1-r) 4 



21(1 - r) 4 + 847(1 - r) 4 j 4 - 16j 3 (l - r) 2 ( - 7(1 - r) 2 + 9 7 (3 + r(-l + 3r))) 

- 4j( - 2(1 - r) 2 (38 + r(-l + 38r)) + 9 7 (l3 + r(8 + r(18 + r(8 + 13r)))) 
+ 12. 7 2 ( - (1 - r) 2 (31 + 7 (-22 + 31r)) + 4 7 (l7 + 7 (-8 + r(27 + r(-8 + 17r)))))] } , 
where in the first step we used that fact that 

~w 4 {j,s)((2s + iy ' 



Y{s) 



(B8) 



(B9) 



APPENDIX C: CONTRIBUTION FROM THE AXISYMMETRIC MODES 

The n = mode, A„ l0 , can be treated independently as follows. Using the binomial expansion, we straightforwardly 
obtain 

(2 5 ) 2(1 



2tt( S - 1) ^ 

(2g)^s) ~ 2-^T(i + 5-l) 



J=0 



1 /l\2(l-s-j) 



= J dx ^^2^ -IfM [C*(2- + 5fc-2,-)-(-J 



(CI) 



This expression requires no further rcgularization, because the analytic continuation is already built-in to this type of 
zcta function. The final result is 



V(», = / 



Cfl(3) 



where we have used the fact that 



E 

j=2 



2 -2 j 

i(i - 1) 



Ch(~2 + 2j, \) - C-) 2 ^} = \{ - 1 + ln2) + ^C*(3) , 



(C2) 



(C3) 
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and Cr(-2) = -C«(3)/47r 2 , which is derived from £ R (z) = it*' 1 / 2 ^! - z)[T((l - z)/2)/T(z/2)]. An important 
observation is that (C2) does not agree with the corresponding expression obtained via the WKB approximation, by 
a factor of —3/8 (see equation (D36) in Ref. [l[). 
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